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1978 — First 2D simulations by Leboeuf et al. 
Early 80s — First 3D simulations by Brecht, Lyon, 
Wu and Ogino 


Late 80s — Model refinements including FACs, 
ionospheres, higher resolution 


90s — ISTP integrates theory and modeling with 
spacecraft missions and comparisons with in situ 
space and ground observations begun 

Today — Global modeling has become integrated 
part of many experimental studies and we’ve 
begun linking models of different regions together 


Ideal MHD Equations 


Non Conservative Formulation 
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e No strict numerical conservation of energy and momentum 


e Various numerical issues 
— Errors in propagating strong shocks 
— Errors in RH Conditions 
— Incorrect shock speeds 


Ideal MHD Equations 


Full Consesyatve Formulation 
oP y pu |=0 


d (pu), 


ot le 
a BP e a 
OO vipero j u. B |- 
87 


V.B =0 


| 1 
=(y-1)| E-— pu” -— B° 
p=(y ) zP" 一 | 


e Strict numerical conservation of mass, momentum and 
energy 
e Numerical difficulties in low (¥ 


— negative pressures possible because p becomes difference of two 
large numbers 


Ideal MHD Equations 


Gas Conservative Formulation 
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。 Strict numerical conservation of mass, momentum and 
plasma energy 
— Final term in Energy equation is j> E dotted with velocity 
— no strict conservation of total energy 


e No difficulties in low (¥j regions 


Time Differencing 


Explicit time differences 
— Predictor — Corrector (2nd order accurate) 
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— Leap Frog Scheme (2nd order accurate) 
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Stability Criterion (CFL Condition) 
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Implicit Schemes generally not used because 
soloution of large linear systems becomes too 
expensive 


MHD Numerics 


e Need a method to solve the conservative 
formulation of the MHD equations which 
maintains the conservation properties 


e For this disucssion we’ ll consider the linear 
advcetion equation 
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Spatial Discretization 


。 State variables are cell 
centered quantities and we 
discretize our model 
equation with numerical 
fluxes through the cell 
interfaces 
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Scheme is conservative 
because 
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Conservative Finite Difference Scheme 


Donor Cell 


A simple first order 
algorithm 
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e Maintains monotonic 
solution 


e Linear advection problem 
clearly shows diffusive 
character 


Donor Cell 


A simple first order 
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e Maintains monotonic 
solution 


e Linear advection problem 
clearly shows diffusive 
character 


Second Order 


2nd Order 


e A simple second order 


algorithm 
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。 Does not maintain 


monotonic solution 


e Introduces dispersion 
errors as seen in linear 
advection example 


Second Order 


e A simple second order 
algorithm 
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。 Does not maintain 
monotonic solution 
e Introduces dispersion 


errors as seen in linear 
advection example 


。 Combines low order and 
high order fluxes 
n+l n vAt 
u, =u, t (Eaz -F an) 


Fan = > (u; 十 2 )+ sign (A Tu; ye 


n n n n 
* max (0, u.,,—u,|—Bs,|u; — u; ) 


S, = = sign (ul, -u, + sign(u! -W ) 
Limiter to keeps solution 
monotonic 
Provides nonlinear 
numeric resistivity and 
viscosity 


Partial Interface Method 


。 Combines low order and 
high order fluxes 
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Limiter to keeps solution 
monotonic 
Provides nonlinear 
numeric resistivity and 
viscosity 


Treatment of the Magnetic Field 


。 Various approaches can be used to satisify the 
constraint that (x)| 


Pi 


— Projection method 
Vv =-VeB 
'=B+VY¥ 


jB convection 
e Modify the MHD equations so that [¥][¥]B convects through 
the system 
d(V eB) 
ee 
dt 


— Use a magnetic flux conservative scheme that keeps 


Magnetic Flux Conservative 


Scheme 
e Magnetic field placed on 
center of cell faces 


Electric field 1s placed at 
center of cell edges so that 


e Cancellation occurs when 
field components of all six 
faces are summed up 


Computation Grids 


。 Simulation boundaries should be in 
supermagnetosonic flow regimes 
— > 18 Re from Earth on Sunward side 
— > 200 Re in tailward direction 
— > 50 Re in transverse directions 


e A variety of grid types exist with varying 
degrees of complexity 
— Uniformed Cartesian 
Stretched Cartesian 
Nested Cartesian 
Regular Noncartesian 
Irregular Noncartesian 
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e Regular Noncartesian 


Domain Decomposition 
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ghost cells 


e Computational Space is divided (evenly?) 
amongst the CPUs available to work on the 
problem 

— MPI used to pass boundary information between ghost 
cells at interfaces 
— Can also use packages like MultiBlockParti and P++ 


Boundary Conditions 


e Upstream 
— Fixed or time dependent values for 8 plasma 
parameters 


。 Can be idealized for derived from solar wind 
observations 


— Problem with B, 
e Need to know 3D structure of solar wind because 


VeB=05n0(B, cam B 


stream downstream 


)=0 
e Implies B,=B, cannot change if solar parameters 
are independent of Y and Z 


e Find n direction with no variation and then sweep 
these fronts across front boundary 


Boundary Conditions Il 
。 All other sides 


— Free flow conditions for plasma and transverse 
components of B 


oF 9 


on 
— normal component of B flows from {Xj x] 
e Inner Boundary Condition 
— MI Coupling module 


— Hard wall boundary condition for normal 
component of velocity and density 


Inner boundary of MHD domain is placed between 2-4 
R, from the Earth 


Parameters in MHD region are mapped along static dipole field 
lines into the 1onosphere 
Field aligned currents (FACs) and precipitation parameters are 
used to solve for ionospheric potential which is mapped back to 
inner boundary as boundary condition for flow 
(-V®)xB 
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Ionosphere Model 


。 2D Electrostatic Model 


îi =0 at low latitude e boundary of ionosphere 
。 Conductivity Models 


—~ Solar EUV ionization 


Creates day/night and winter/summer asymmetries 
2, =0.5Fo; (cosy) V x <65° 
>, =1.8F 4; cos y V ps<65% 


- Auroral Precipitation 


Empirical determination of energetic electron precipitation 


Auroral Precipitation Model 


e Empirical relationships are used to convert MHD 
parameters into an average energy and flux of the 
precipitating electrons 


— Initial flux and energy 
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— Parallel Potential drops (Knight relationship) 
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— Hall and Pederson Conductance from electron precp (Hardy) 
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e A whole lot of coding later and you get 


Methods of Model Validation 


Conduct studies with same conditions and different 
numerics 


Computation of theoretical problems with known 
analytic answers 
— Provides a ground truth that code is working 


— Very limited number of MHD problems 


Direct comparison with observations 


— Limited number of spacecraft observations 
e Check general characteristics with superposed epoch studies 
— Include comparison with indirect observations 


— Use metrics to quantitatively asses validity 
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。 Low Numerical Diffusion 


e High Numerical Diffusion 
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e Simulation for Northward IMF with constant Pedersen Conductance 


e High Numerical Diffusion 。 Low Numerical Diffusion 


。 Both data and Simulation show onset, 
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intensification, recovery, and second onset 


。 Simulated onset is early, but intervals between 


intensification and second onset are consistent 
Simulated CL recovers faster than observations 


Medium Skill Case 1 Raeder: 0.27 


DMSP data-model comparison 。 ; 
a1398310.154 Ridley: 0.32 


Rice: 0.32 
Slinker: 0.23 
Weimer: 0.43 
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Spatial structure details not reproduced 
More detail (Rice, Slinker) does not necessarily lead to higher skill, 
although it is desirable 


Comparison with geostationary 
observations 


。 Excellent agreement for 


all three components of 
B 


Despite global B, 

offset dipolarizations of 

similar size are seen in 

simulation results for 

both GOES 8 & 9 

— May imply limited role 

for ring current in 
substorms 


Flow Channels 


Comparison between 
Flow channels and BBFs 


INNER CPS FLOWS 


e Flow channels have 
properties similar to BBF 
results reported by 
Angelopolous 
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FWHM of V,, profile and 
magnitude comparable 
BBF properties 


Use code to determine if 
they result from localized 
reconnection or 
interchange instability 


LFM-TING Coupling 


Coordinate transfer 
Data interpolation 


Particle precipitation: F,, Eo 


One Way Coupling 


Conductances 
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Two Way Coupling 


Conclusions 


。 Global MHD simulation of the magnetosphere under 
idealized solar wind conditions are proving to be a 
useful tool for expanding our understanding the 
coupled solar wind — magnetosphere — ionosphere 
system 


e The technique is expanding into new frontiers 


— Jonospheric simulation is being replaced with more 
sophisticated Thermosphere-Ionosphere Global Circulation 
Models 


— Modeling of the inner magnetosphere is being enhanced by 
coupling with the Rice Convection Model 


Conclusions 


e LFM is highly successful global MHD simulation of 
the magnetosphere 
— Numerous publications and presentations 
— Its design considerations are still relevant today 
e LFM is still evolving 
— Ports to new platforms and utilization of MPI 


— Jonospheric simulation is being replaced with more 
sophisticated Thermosphere-Ionosphere Global Circulation 
Models 


— Modeling of the inner magnetosphere 1s being enhanced by 
coupling with the Rice Convection Model 


